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ABSTRACT 

This paper provides full sky maps of foreground emission in all WMAP channels, 
with very low residual contamination from the Cosmic Microwave Background (CMB) 
anisotropics and controlled level of instrumental noise. Foreground maps are obtained 
by subtraction of a properly filtered CMB map, obtained from linear combinations of 
needlet-based representations of all WMAP observations and of a 100-micron map. 
The error in the reconstructed foreground maps on large scales is significantly lower 
than the original error due to CMB contamination, while remaining of the order of 
the original WMAP noise on small scales. The level of the noise is estimated, which 
permits to implement local filters for maximising the local signal to noise ratio. An 
example of such filtering, which reduces the small scale noise using latitude dependent 
filters is implemented. This enhances significantly the contrast of galactic emission, in 
particular on intermediate angular scales and at intermediate galactic latitude. The 
clean WMAP foreground maps can be used to study the galactic interstellar medium, 
in particular for the highest frequency channels for which the proper subtraction of 
CMB contamination is mandatory. The foregrounds maps can be downloaded from a 
dedicated web site. 
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1 INTRODUCTION 

The WMAP space mission, launched by NASA in 2001, has 
been primarily designed to measure the anisotropics of the 
Cosmic Microwave Background (CMB) emitted when the 
universe became transparent, at an age of about 380,000 
yr. The interpretation of these observations in the context 
of the standard hot Big-Bang model has constrained the 
main parameters of the model with great accuracy (Komatsu 
et al., 2009). 

In addition to the primordial CMB, WMAP detec- 
tors are sensitive to foreground astrophysical emission: dif- 
fuse emission from the galactic interstellar medium (ISM), 
emission from compact extragalactic sources (such as radio- 
galaxies and AGN and, to lesser extent, thermal Sunyaev- 
Zel'dovich (SZ) emission from clusters of galaxies). Hence, 
maps observed by WMAP contain each a mixture of emis- 
sions from different astrophysical processes. 

The analysis of the observed maps for deriving a CMB 
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power spectrum, in temperature or polarization, involves 
some kind of foreground cleaning, followed by masking the 
regions of the sky most contaminated by foregrounds, and 
estimation of the angular spectrum Cf^^ of the CMB on 
the masked sky. There exist several ways to deal with the in- 
complete sky coverage for power spectrum estimation (Hivon 
et al., 2002; Efstathiou, 2004; Fay et al., 2008). Little is lost 
in terms of power spectrum estimation accuracy by masking 
a small fraction of sky (the accuracy is typically reduced by 
a factor /sky equal to the fraction of the sky kept for the 
analysis). Proper estimation of the contamination (biasing) 
of the power spectrum by noise or residual foreground emis- 
sion permits to correct for the contribution of both kinds of 
contaminants on the power spectrum of the CMB map. 

Where to give up foreground cleaning and resort, in- 
stead, to masking, is an interesting (and debated) issue. To 
be on the safe side, conservative masking is typically used 
for estimating the CMB power spectrum. However, low-level 
foreground emission often remains in the masked sky (see, 
e.g., Patanchon et al. (2005); Leach et al. (2008)). 

Foreground emissions of various astrophysical origins, 
however, are of much scientific interest themselves. Scien- 
tific investigation of foregrounds cannot resort to masking 
regions contaminated by CMB, the latter being present ev- 
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ery where on the sky. In recent observations, such as those 
of WMAP, the CMB is stronger than the noise on a large 
range of angular scales, in all regions of the sky. The analysis 
of foreground emission in WMAP data, hence, requires sub- 
traction of the CMB contaminant (which is the main source 
of error) while keeping contamination by instrumental noise 
as small as possible. Obtaining such clean foreground inten- 
sity maps from WMAP data is the objective of the present 
paper. 



2 METHOD 

Standard methods for component separation (see De- 
labrouille & Cardoso (2007) for a review) often assume that 
the data are well represented by noisy linear mixtures of 
well-defined components, e.g. CMB, SZ, thermal dust, syn- 
chrotron, etc. Most methods used in the CMB context such 
as Wiener filtering (Bouchet & Gispert, 1999; Tegmark & Ef- 
stathiou, 1996), Maximum Entropy Methods (Hobson et al., 
1998), Independent Component Analysis (ICA) methods of 
different types (Maino et al., 2002; Delabrouille et al., 2003) 
explicitly represent sky emission as a superposition of emis- 
sions, the properties of which are parametrised in some way. 
In reality, matter emits via a great variety of distinct pro- 
cesses, and there is no obvious way of separating the total 
emission into distinct components. It is probably fair to ar- 
gue that there is, for instance, no natural preferred option 
between distinguishing components via their emission pro- 
cess (e.g. synchrotron vs. free-free) or via their place of ori- 
gin (the ISM of our galaxy vs. the ISM of other galaxies). 
Similarly, dust emission arises from a variety of forms of 
matter (from molecules to large dust grains), and from sev- 
eral emission processes (thermal greybody emission, electric 
dipole emission from rotating dust grains...). Whether to 
distinguish thermal dust from spinning dust as two separate 
components, or instead distinguish between warm and cold 
dust, for instance, is a matter of taste. One can also consider 
dust emission as being all emission from all matter that is 
not ionised - the appropriate modeling depends on the sci- 
entific question addressed. The obvious exception to this is 
the CMB, which is distinct both by its origin (the last scat- 
tering surface) and its spectral emission law. If, however, the 
WMAP experiment was sensitive enough for kinetic SZ ef- 
fect (which has the same spectral emission law as primordial 
CMB fluctuations) to be above the noise level, the question 
of whether to distinguish the two would become relevant as 
well. 

2.1 The model 

Given that complexity of foreground emission, and the rela- 
tive lack of spectral resolution in WMAP data (only 5 chan- 
nels, for many different emission processes), one must extend 
the way to parametrise foreground emission. Extensions of 
the Spectral Matching ICA (SMICA) method (Delabrouille 
et al., 2003) to allow more flexible models of foreground 
emission are discussed by Cardoso et al. (2008a, b). 

In the present analysis, we adopt the extreme point of 
view in which we choose to distinguish between three contri- 
butions to the WMAP maps only: CMB, foregrounds (dom- 
inated by the emission of the galactic ISM), and noise. The 



task we address, then, is to produce the best possible maps 
of foregrounds with minimal contamination from CMB and 
noise. Further analysis of the foreground component into dif- 
ferent emissions of interest (whether by origin, or by emis- 
sion process) is not the aim of this paper. 

The temperature map at a given frequency channel of 
WMAP is given by 

x{p) = j h\p.p') [s{p') + r (pO] dp' + n\p), (1) 

where i indexes the WMAP frequency bands, and runs from 
1 to 5. s{p) and /*(p) are respectively the CMB and fore- 
ground components, and n*(p) is the instrumental noise in 
the given frequency band. Here, b^{p.p) is the beam of the 
observed map for the i^^ frequency band (assumed to be 
symmetric, but not necessarily Gaussian). The direction on 
the sky (or pixel of a pixelized map) is indexed by p. Equa- 
tion 1 can be written in harmonic space as: 

Xlm ^Sim + /Im] + ^^m (2) 

Our objective is to estimate in the best possible way the fore- 
ground emission f b^{p-p)f^{p)dp (or equivalently h\fl^) 
in each WMAP channel. Once stated in that way, the prob- 
lem becomes quite well posed. The CMB emission has the 
particular feature that it is completely uncorrelated from 
the other two, and that its emission law is known (it is the 
derivative of a blackbody with respect to temperature, at 
temperature Tcmb = 2.725 K). The instrumental noise is 
well characterised. It is well approximated as uncorrelated 
from channel to channel. Its level is known with excellent 
accuracy, and it is known to dominate at high angular fre- 
quencies (high I). 

The foreground component is everything else, i.e. ev- 
erything that comes from the sky (rather than from the 
instrument) and is not CMB. We set no other constraint 
on this foreground component, i.e. we assume no particular 
parametrization of it. 



2.2 Separation strategy 

We propose to extract the foreground component by sub- 
tracting an estimate of the CMB from the total maps, and 
then post-processing the remaining set of CMB-cleaned 
WMAP data to get maps of foreground emission. The ef- 
fectiveness of this procedure is discussed in section 4 and 
appendix A. 

Considering the minimal assumptions made about the 
foregrounds, it is quite natural to use the so-called Inter- 
nal Linear Combination method (ILC) to extract a CMB 
map, to be subtracted from WMAP data to produce CMB- 
free foreground maps. The ILC is a multifrequency linear 
filter which uses a linear combination of the input data. It 
minimizes the variance of the reconstructed CMB map over 
domains of observation as, e.g, regions of the sky, or domains 
in harmonic space, or domains in needlet space. 

The strategy for obtaining foreground maps is thus the 
following: 

• estimate a CMB map using an ILC; 

• subtract that CMB from the WMAP frequency maps 
(at proper resolution); 
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Figure 1. The clean CMB map extracted using internal linear 
combination in needlet space, used for CMB subtraction in the 
present analysis 



• filter to maximize the signal to noise ratio in each of 
the foreground maps. 

2.3 The needlet ILC 

Considering the variability of the statistical properties of 
foregrounds and noise both in pixel space and in harmonic 
space, the optimisation of the separation calls for locali- 
sation of the filter in both spaces. That can be achieved 
with needlets, which are a special type of wavelets on the 
sphere (Narcowich et al., 2006; Marinucci et al., 2008; Guil- 
loux et al., 2009). A needlet-based ILC is well suited to 
our component separation problem, as localisation in direct 
space (in addition to localisation in harmonic space) per- 
mits better reconstruction of the CMB in the vicinity of the 
galactic plane, i.e. a region of major interest for foreground 
science. For the present analysis, we use for CMB subtrac- 
tion the WMAP 5-year needlet-based ILC map (NILC map) 
obtained by Delabrouille et al. (2009), and displayed in fig- 
ure 1. Other options were considered, as for instance the 
harmonic ILC map derived in Saha et al. (2006, 2008), as 
has been done in Ghosh et al. (2009). The NILC map has 
been chosen here because it reconstructs best the CMB in 
regions where foregrounds are significant. It can be written 
as: 

S£m — Sim + ^im (3) 

where 6im is the CMB reconstruction error and where the 
beam is denoted bf^ , corresponding to the resolution of chan- 
nel W. The angular spectrum of that map is 

E|s,^|^ = Brcf^^+iVf^ (4) 
with A^iLC ^ E{Sem)'' and = (b^^f. 

2.4 Channel— dependent CMB subtraction 

In our present application, the objective is to subtract the 
CMB in each channel, without adding significant additional 
noise (and, obviously, without adding more noise than we 
subtract CMB). In each channel, the total emission observed 
(eq. 2) comprises a CMB contribution blsim- We subtract 
from each observation map x}^ a rebeamed, rescaled version 



of the needlet ILC map, that is, the foreground map fl^ at 
the resolution of channel i is estimated by 

hi 

/i i '^i ^ /r\ 

where the coefficient ai can be chosen to minimise the re- 
construction error, denoted s}^, and defined by: 

fim — blfim -\- £im- (6) 

Indeed, eqs (2), (3) and (5) show that 

b^ 

^Im = (1 - Oii)blSim - ai-^Sim + Tl^rn (7) 

Neglecting the correlation between 6irn and n^^, one gets 

E(££m) =(!-«€) +^^-^^1 +E(n^^) (8) 

where Bl — {b\Y . The mean square error E(££^)^ is easily 
found to be minimal for 

^CMB 

This is exactly the signal to signal-plus- noise ratio in (4), 
that is, the prescription for Wiener filtering. Therefore, the 
variance of the foreground map in channel i (at the resolu- 
tion of the channel), is minimal if we subtract the Wiener- 
filtered version of the needlet ILC map. 

CMB subtraction according to (5) yields our primary 
foreground map in WMAP channel i, at the resolution of 
the considered WMAP channel. The map contains essen- 
tially foreground emission and noise from channel i, and 
additional error originating from imperfect CMB subtrac- 
tion. This additional error is discussed in more detail in sec- 
tion 4.3. Figure 2 illustrates for one of the WMAP channels 
(here, the K band at 23 GHz), on a patch centered at modest 
galactic latitude (latitude -30 degree and longitude 70 de- 
gree), the separation of the original observation into CMB 
and foregrounds. The effect of CMB subtraction is clearly 
visible when comparing the original map (top left) to the 
map after CMB has been subtracted (bottom left). Most of 
the large to intermediate scale fiuctuations due to the pres- 
ence of CMB have been removed, which permits to identify 
clearly compact sources and structures of the galactic emis- 
sion. 

2.5 Latitude-dependent filtering 

After CMB subtraction, the estimated foreground map fijji 
still contains significant noise, which dominates on small 
scales. This is clearly visible in the bottom left panel of figure 
2. That map is obtained, by construction, at the resolution 
of the original input map but we still have extra degrees 
of freedom: adjusting the beam or, equivalent ly, multiplying 
the harmonic modes by some filter one obtains a new 
estimate f}^ of the foreground map as: 

fi i pi 

Jim — U)£j£m 

The filter can be adjusted to minimize the total variance 

^ , . . . ,2 

^Ifim—^ifiml 5 
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and it is well known that, in the ideal case where both the 
foreground map of interest and the contamination term 
are stationary, the best filter is the Wiener filter: 

~ BICJ^ + Et 

where Ei is the variance of the reconstruction error e\^. 

However, we are not facing a stationary situation. As 
galactic foregrounds are strongly concentrated in the galac- 
tic plane, the filtering should better depend on the location 
on the sky. It is not immediately clear what would be the 
best strategy for such a localized filtering, especially if there 
is need not only to provide the best map but also to char- 
acterize it. As a reasonable compromise between efficiency 
and simplicity, we approximate the foreground emission as 
following a plane parallel slab model. We then apply a filter 
for each zone of galactic latitude. As discussed in the next 
section, we also constrain the total response after filtering 
to have a Gaussian shape, i.e. we find the Gaussian beam 
which fits best the harmonic response w\h\ resulting from 
the instrumental beam and the Wiener filter. In that way, 
the variation of the beam across the sky can be described 
simply by a latitude dependent (actually: depending only 
on the latitude zone) Gaussian beam. The improvement of 



the quality of the derived foreground map after filtering is 
clearly visible in figure 2, where the bottom right panel dis- 
plays the foreground map after the filtering. Details of the 
implementation are given in the next section. 



3 PRACTICAL IMPLEMENTATION 
3.1 CMB subtraction 

For each WMAP channel, we start from the Wiener filtered 
Needlet ILC CMB from Delabrouille et al. (2009) which in- 
cludes the ai ratio of eq. 9. The resolution of the map is 
changed in harmonic space, by multiplication of each a^m 
by the ratio h^i/hf . The resulting CMB map is subtracted 
from the observation in channel i. Note that because of the 
Wiener filter applied on the CMB map, the equivalent res- 
olution of this CMB map is not that of channel i. 

The foreground map obtained in a given channel by 
subtraction of the CMB contains foregrounds, noise, and 
a residual (the difference between the CMB estimate and 
the true CMB at the resolution of the given channel). The 
residual contains some noise from all channels, some resid- 
ual CMB, and a small amount of foregrounds. This resid- 
ual, however, is well below the foreground emission on large 
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Figure 3. Angular power spectra over the full sky for the frequency bands K to V. The smooth black curve is the power spectrum of the 
WMAP CMB-subtracted map, the red line (dashed) is the theoretical WMAP noise spectrum for the channel considered, and the orange 
line (double dashed) is the noise contribution from the NILC map. Foreground (galactic) emission is dominant on large scale, whereas the 
theoretical WMAP noise nearly matches the power spectrum of the reconstructed map on small scale. The NILC noise lies somewhere 
between the level of the CMB and that of the WMAP channel noise on large scales, and is significantly below the instrumental noise on 
small scales. 



scales and orders of magnitude below the noise of the con- 
sidered WMAP band on small scales. This is illustrated, 
for four of the WMAP channels, in figure 3, which displays 
the power spectra of the WMAP observations after CMB 
subtraction, as compared to the original CMB contamina- 
tion, to contamination by noise, and to the residual due to 
improper CMB subtraction. In all channels, the process of 
CMB subtraction reduces by two orders of magnitude the 
total contamination power on the largest scales. For the V 
channel, a significant reduction of the total contamination 
is achieved up to i larger than 500. 

3.2 Latitude zones 



Here is the width of transition regions of the mask (re- 
gions in which its value decreases smoothly from 1 to 0) . The 
latitudes Oq and Oi are the lower and upper limits of the re- 
gion where the mask value is 1. Note that each region is 
symmetric with respect to the galactic plane. Smooth tran- 
sitions in the masks for the different latitude bands simply 
permit to avoid edge effects in the masked maps, mitigating 
mode mixing in the calculation of power spectra. Limits for 
the zones (i.e. values for and ^i) used in this analysis 
are 0, 5, 10, 20, 30, 45, 60 and 90 degrees, with transitions 
of = 5°. The width of the regions selected in this way 
increases with increasing galactic latitude, which is required 
for accurate estimation of the Wiener filter (see below). 



The sky is divided into 7 latitude zones, selected using mul- 
tiplicative masks. Each mask is defined as follows: 



M(6',(/)) = 1 



: 



for 6*0 < \e\ < 01 



= cos' f I X ) for ^0 - A^ ^ 1^1 < Oo 

: cos^ X ) for ^1 < 1^1 < ^1 + A^ 



elsewhere. 



3.3 Latitude-dependent Wiener filter 

For each channel, for each latitude band, we compute the 
power spectrum C^^^ of the masked noisy foreground map. 
For each i, the total power comprises a contribution Fi from 
the foregrounds b}/}^ at the resolution of channel i (i.e. the 
signal of interest) and a contribution Ei from the recon- 
struction error e}^. In this case, the Wiener filter can be 
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Figure 4. Angular power spectra after the Wiener filtering for the frequency bands K to V. This is the equivalent of fig. 3 after noise 
filtering. The bump visible on small scales for channel K (and, to lesser extent, channel Ka), is due to the non-gaussian shape of the 
WMAP beams. It could be gotten rid of by relaxing the constraint that our final map should be characterised by a Gaussian beam. 



approximated as: 

^OUT 77. 

-^ = %pw^- (10) 

See section 4.3 for the estimation of Ei. The Wiener filter 
as defined by equation 10 is applied independently in each 
WMAP channel. 

The effective beam after Wiener filtering becomes wib}. 
That is the beam which minimizes the variance of the error 
of the reconstructed foreground map, for channel i, in the 
latitude band defined by the mask corresponding to the zone 
considered. In fact, such a Wiener filter can be computed 
locally for any region of interest, not only for latitude bands, 
which are used here mainly for illustrative purposes. 

3.4 Gaussian beam approximation 

As discussed above, it is not very convenient to construct a 
final foreground map in which the beam varies rapidly over 
the sky. For simplicity, it may also be desirable to work with 
effective Gaussian beams, in which case the beam is fully 
defined by a single number: its full width at half maximum 
(FWHM). The beam response in harmonic space then is: 



where a is related to the FWHM of the Gaussian beam by 
the relation a = ^fwhm/ v^(8 x log 2). 



For each of our seven selected regions then, we deter- 
mine the Gaussian beam which best matches the theoretical 
optimal effective beam due to the effect of both the instru- 
mental beam and the Wiener filter. We derive a FWHM 
which varies smoothly in latitude bands of 5 degree width by 
interpolating the FWHM derived using the broader latitude 
bands described in 3.2 above. The corresponding FWHM, 
for all WMAP channels and all latitudes in 5 degree steps, 
is plotted in figure 5. 



Note that while our choice results from a trade-off be- 
tween optimality and convenience, it is somewhat arbitrary. 
We could, for instance, have decided to work with more com- 
plex symmetrical beam shapes or different zones. We stress, 
however, that no information is lost in that process since 
one can always change back the response in harmonic space 
(and change the effective beam) by multiplication of a£rn 
coefficients by the appropriate response ratio. The FWHM 
plotted in figure 5 are just indicative of the scales where 
noise starts to dominate over the foregrounds, as a function 
of channel and galactic latitude. 
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Figure 5. The FWHM of the Gaussian beam used to approxi- 
mate the Wiener filter, as a function of latitude, for all WMAP 
channels. For each WMAP frequency channel, the expected gen- 
eral behaviour is observed: the resolution can be better at low 
galactic latitude, where there is more signal, and hence better 
signal to noise ratio at high £. 



3.5 Construction of the filtered maps 

After having determined, for each frequency channel, the 
resolution at which the foreground map should be recon- 
structed as a function of galactic latitude, a filtered final 
map for each WMAP channel is obtained in the following 
way. Starting from the original WMAP observations and af- 
ter CMB removal, we make 18 different maps, one at the 
appropriate resolution for a 5-degree wide latitude band. 
The final map for that channel is obtained by picking, for 
each pixel of the final map, the corresponding pixel in the 
one of the 18 maps which is at the resolution matching the 
latitude band of the pixel at hand. Hence, the resolution of 
the filtered map varies in steps over the sky, with a change 
in resolution at all multiples of 5°. 



4 DISCUSSION 

4.1 Mult i- frequency analyses 

As seen in figure 5, the Gaussian beam which maximises the 
signal to noise ratio of the foreground map in a given channel 
depends both on galactic latitude and on the channel. 

Further analysis of these foreground maps, however, in- 
volves measuring spectral indices in selected regions or the 
spectral energy distribution of compact sources, i.e. ratios 
between flux at different frequencies. Such post-analyses re- 
quire maps at the same resolution. 

For speciflc, localised analyses, using several WMAP 
channels, original CMB-cleaned maps can easily be re- 
beamed for better matching of the resolution with the local 
S/N ratio. Figure 5 provides a helpful tool for choosing the 
resolution of the analysis. To take a speciflc example, for in- 
vestigations involving all channels at 25 degrees latitude, one 
should pick a resolution somewhere between 30 and 50 arc- 
minutes, as the curves giving the appropriate resolution, for 
each channel, in that latitude zone, give numbers between 
30 and 50. 



4.2 Faint foregrounds? 

Although the original claim in the present work is the con- 
struction of maps of all foregrounds, it should be noted that 
the flnal maps are appropriate only for the scientiflc analysis 
of foregrounds whose contribution to the observed emission 
is larger than the WMAP noise (but can be smaller than the 
CMB), i.e. of galactic foregrounds and strong point sources. 

The reason is as follows. The subtracted CMB has been 
obtained by an ILC method, which flnds locally linear com- 
binations which minimise the total variance of the CMB 
reconstruction error. When subdominant foregrounds are 
present, it is more favourable to the ILC to let weak fore- 
grounds leak in the CMB map and reduce the instrumental 
noise, than the opposite. For instance, Sunyaev Zel'dovich 
effects are too weak to play a role in the adjustment of the 
linear combinations used for CMB reconstruction in WMAP 
data. The same is true for a weak background of point 
sources. Hence, the ILC is not optimised to reject these fore- 
grounds, which thus are not preserved in the flnal foreground 
maps produced with our processing. A search for SZ signals 
in WMAP data has been performed by Melin et al. (2010) 
using a method appropriate for this particular signal, the 
multifrequency matched fllter (Herranz et al., 2002; Melin 
et al., 2006). 



4.3 Noise levels 

The 'noise' in our foreground maps (from all sources of error) 
in the map for channel i (at the resolution deflned by the 
beam b}) comes from four main terms: 

(i) noise from the relevant WMAP channel; 

(ii) noise "^wju^ from 'internal linear combinations' of 
WMAP channels, present in the CMB map subtracted from 
each WMAP observation; 

(iii) foregrounds leaking in the needlet ILC estimate of 
the CMB; 

(iv) CMB residuals after subtraction of the estimated 
CMB. 

An exact characterisation of the total error is not completely 
possible in the absence of a reliable model of the foregrounds, 
which is needed to estimate the amount of foreground emis- 
sion leaking in the estimated CMB. It is possible, however, 
to give an estimate of the total 'noise' power spectrum, E^. 



4.3.1 Noise power spectrum estimate 

An estimate is obtained by considering that all terms are 
independent so that the total power of the error is the sum 
of the powers of all contributions. 

The flrst term is simply the noise power of the 
WMAP i^^ channel, the level of which is provided with the 
release of WMAP data. 

The second term is estimated by Monte- Carlo simula- 
tions of the WMAP noise by Delabrouille et al. (2009). Noise 
maps and power spectra released are available on a dedicated 
web page^. They should be corrected from the effect of the 



http://'w 



-paris7.fr/ APC/ Recherche /Ad amis /cmb_wmap-en.php 



8 T. Ghosh et. al 




-&0,0) Galoctio 



Figure 6. 40 degree by 25 degree patch of the sky at latitude of -50 and longitude of +200. The K and Ka band foreground maps are 
shown in this range. Most of the point sources (circled by black) identified by the WMAP team are visible by eye. 



Wiener filtering and rebeaming. This is done by multiplying 
each mode, in harmonic space, by: 



Ka Band 



(11) 



where wi is the effective beam of the Wiener-filtered needlet 
ILC CMB map, b} is the beam of the foreground map, and 
b^^^^ the beam of the NILC map (coinciding, for i < 1000, 
with that of the W channel. The power iV,c^^ of the error 
from this term is: 

iV^ ^riNi (12) 

where A/"]^^ is the noise level, per harmonic mode, available 
with the needlet ILC CMB map. 

The third term is unknown. On the vast majority of the 
sky where foregrounds are weak, it is negligible. Where fore- 
grounds are strong, it is always a small fraction of the input 
foregrounds. It is therefore neglected in our noise estimate. 

The fourth term comprises contributions from two ef- 
fects: a) residual CMB from the difference between the beam 
of the WMAP channel considered and the effective beam of 
the subtracted CMB, which can be computed straightfor- 
wardly; b) CMB modes cancelled by the ILC because of 
empirical correlation with the foregrounds and noise (see 
appendix and figure 6 and of Delabrouille et al. (2009)). 

Appendix B gives a detailed calculation of all noise 
terms. The power of the sum of all contributions (to- 
tal error) is given by the expression: 



(bl (1 - ai)y + OMbl (1 - ai) bfn 



C?^^ (13) 



where C^^^ is the CMB power spectrum (which can be 
taken to be the WMAP best fit model, for instance). The 
last term is the total contribution from CMB subtraction 
error. The first term in the brackets comes from effect a) 
above, and the second term is a correction originating from 
the 2% CMB power loss due to the ILC bias. 

Figure 7 shows the decomposition of the total noise into 
three main terms: noise from WMAP in channel Ka, noise 




Multipole(0 

Figure 7. Power spectra of the noise in the Ka band, decomposed 
into three main sources of error. The instrumental noise from 
that frequency channel dominates on small scales, whereas CMB 
reconstruction error (noise in the NILC map) dominates on large 
scales. The contamination of the foreground map by the CMB is 
very small. 



from the NILC map, and residual CMB improperly sub- 
tracted. 

Figure 8 displays the power spectra of all foreground 
maps, and power spectra of the noise term in all of them. 



4-3.2 Noise power spectrum upper limit 

The two first terms (i) and (ii) (arising from WMAP instru- 
mental noise) are correlated because the instrumental noise 
in the estimated CMB contains a part of the noise on the 
channel of interest. An upper limit to the total error is ob- 
tained by assuming that the correlation is perfect between 
and "^wju-^ , which is the case when channel i dominates 
with a negative coefficient in the ILC (negative, because the 
ILC is subtracted from channel i to get the foreground map). 
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Figure 8. Power spectra of foreground maps and corresponding 
noise level. 



Hence, we can assume the following upper limit for Ei 

2 



+ 



(bl (1 - a,))' + 0.04 6^ (1 - ai) b^n 



C?^^ (14) 



4-3.3 Noise level maps 

The origin of all main contributions to noise being identi- 
fied, it is a straightforward matter to write a pixel space 
equivalent of equations 13 and 14. 



4.3.4 More accurate noise characterisation ? 

The description of noise as a single noise power spectrum 
per map only is an approximation. It is convenient for most 
applications, but one should be aware of its limitations. 

The first limitation comes from the fact that the noise 
is not stationary for two reasons. The first reason is that 
the WMAP hit count is not homogeneous on the sky. The 
second reason is that galactic foregrounds, and hence the co- 
efficients of the needlet space ILC used to subtract the CMB, 
and hence the noise from the ILC map, are non stationary. 
If needed, it is possible to scale the two first terms of equa- 
tion 13 as a function of the local noise levels. One can also 
generalise equations 13 and 14 using the pixel-based noise 
variance of WMAP channel i in place of N} , and using noise 
simulations provided by Delabrouille et al. (2009), filtered 
by in place of Nj^^ . 

The second limitation comes from the cross-channel 
correlation of the resulting noise. All terms in the right 
hand side of equation 13, except the first term, give rise 
to such correlated noise. On large scales, these correlated 
terms dominate the total error, so the noise is close to 100% 
correlated. On small scales, the noise from WMAP i'^ chan- 
nel dominates, so that the noise is not correlated between 
channels. 

More accurate noise characterisation are possible, but 
they require Monte-Carlo simulations of the whole process. 



If needed for very precise analyses, such simulations (in- 
cluding in particular cross correlation of the error between 
channels) can be performed by the authors upon request. 



4.4 Products 

Our foreground maps and related data comprise one CMB- 
cleaned map (i.e. a foreground map) for each WMAP chan- 
nel, at the original resolution of the channel, characterised 
by a corresponding noise power spectrum (given by equation 
13) for each channel and an effective beam for each chan- 
nel (copied from the beam provided with the WMAP 5-year 
data). 

For each channel, a suggested resolution (Gaussian full 
width at half maximum) as a function of galactic latitude, 
appropriate for minimizing the foreground reconstruction er- 
ror at that latitude, is given in figure 5. 

Figure 9 displays our filtered foreground map for chan- 
nels Ka to W in the region of the Gum nebula. Figures Bl to 
B5 display full sky foreground products, both at the original 
WMAP resolution, and with noise filtered out (latitude de- 
pendent beams). Contribution from different astrophysical 
processes can be seen even in the high frequency channels 
on the filtered maps. 



4.5 Comments about the present approach 

The approach to foreground estimation discussed in this pa- 
per is quite possibly the simplest way to achieve reliable 
foreground estimation. 

The foreground products delivered by the procedures 
described in this paper are obtained with little prior infor- 
mation: the input maps are supposed to be well calibrated 
and they are supposed to be accurately characterized by 
their beams and noise properties. Besides that, no prior in- 
formation is used about the foreground emissions, and no 
attempt is made at constraining or modeling those emis- 
sions. As a result, we obtain 'maximal' foregrounds maps. 

It would be interesting to aim at producing more con- 
strained foreground maps by including prior information 
about foreground emission. In the SMICA approach, for in- 
stance, the foreground emission is typically modeled as the 
superposition of a number A^fg correlated templates. This, 
however, is far from being as straightforward as the ILC- 
subtraction proposed in this paper. 

It should be noted that the ILC-subtraction approach 
based on A^obs observation frequencies can be shown to be 
equivalent to assuming that the foreground emission can be 
represented as the superposition of A^fg = A^obs — 1 tem- 
plates (see appendix for a proof and a more careful state- 
ment). Along the same lines, one may also consider a multi- 
dimensional version of the ILC targeted directly at the fore- 
grounds. This is briefly discussed in the appendix and is the 
topic of a forthcoming publication Remazeilles et al. (2010). 



4.6 Comparison with other WMAP foreground 
maps 

Other maps of WMAP galactic foregrounds have been ob- 
tained by a variety of methods. All such maps, however, are 
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Ka-Band Q-Band 




Figure 9. Gum Nebula as viewed from Ka to W band at K-band 
resolution. K band resolution is used to show the diffuse emission 
rather than the small scale features. 
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produced at resolution of 1 degree or worse. The maps ob- 
tained in the present work are provided at the resolution of 
the original WMAP channels. 

In addition, as the method subtracts a GMB map sig- 
nificantly cleaner, in the galactic plane, than that obtained 
by a simple ILG, our maps give better estimates of the total 
emission of the ISM in the vicinity of the galactic plane. 



5 CONCLUSION 

In the present paper, we have analysed the WMAP 5-year 
temperature maps to clean them from the emission orig- 
inating from the cosmic microwave background. The maps 
obtained are noisy foreground maps, where the astrophysical 
emission is dominated by emission from the galactic inter- 
stellar medium, and a small number of compact sources. 

The noise in these maps can be reduced by filtering. 
We estimate, as a function of galactic latitude, the Gaussian 
beam to be used for maximising the contrast of the signal 
of interest in the map. 

The maps produced should not be used for searching 
for faint components such as the Sunyaev Zel'dovich effect 
or the emission from a background of faint sources. They 
are adequate for studying the emission of the galactic in- 
terstellar medium and the spectral energy distribution of 
strong compact sources. Products can be downloaded from 
a dedicated web page.^ 



http:/ /www. ape. univ-paris7.fr / APC/ Recherche / Adamis/fg_wmap-en.php 



http: //lambda, gsfc.nasa.gov/ 
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APPENDIX A: ILC AND ILC SUBTRACTION AS OBLIQUE PROJECTIONS 

In this appendix, we exhibit the direct connection between two strategies for estimating the galactic emission: i) the strategy 
developed in this paper where the ILC-based estimate is subtracted from the observations and ii) a direct estimation of the 
galactic emission based on a multi-dimensional ILC. 

That issue is discussed in the following context. Since the processing takes place in harmonic space, the discussion can be 
simplified by focusing on a single (^, m) mode. Given then A^obs input channels, denote x the A^obs x 1 vector of the spherical 
harmonic coefficients of all channels after beam correction at a particular (i,m). We can write x — as -\- f where s is the 
spherical harmonic coefficient of the CMB for that mode and the A^obs x 1 vector a contains channel gain with respect to 
CMB (it is a vector of I's for a perfectly calibrated instrument and input maps in the appropriate units). Vector / represents 
all the other emissions. 

Al Component estimates as projections 

The ILC estimate of the CMB signal s is s — w^x where the A^obs x 1 vector w is such that w^x has minimum variance and 
has unit gain towards the CMB. Hence, the ILC filter is the minimizer of E{w^xY subject to to* a = 1. It is easily found that 
the solution is — {a^ a)~^ o!' where R denotes the covariance matrix of x. Hence, the CMB reconstruction on all 
channels is the vector as — Pix where matrix Pi is defined by 



am-^a 

Matrix Pi is a projection matrix (Pf = Pi) but it is not an orthogonal projection (Pi ^ Pi). It is an oblique projection 
onto Span(a) along its null space A/" which, by definition, is the (A^obs — l)-dimensional subspace of 

J\r= {z\ a'R-^z = 0} . 

This is the subspace of all directions which are nulled out by ILC. 

In this paper, we considered reconstructing the foreground emission by subtraction: f = x — as so that the foreground 
estimate is related to the data by / = P2X with 

P2=/n-Pl. 

By construction, P2 is the projection matrix onto JV along Span(a). 

Hence, the decomposition of the observed vector into CMB and foregrounds corresponds to oblique projections onto a 
pair (Span(a), A/") of complementary but non orthogonal subspaces of dimensions 1 and A^obs — 1 respectively. 

A2 Multi-dimensional ILC 

One may consider a direct estimation of the foregrounds by generalizing the ILC method to address the case of a multi- 
dimensional component. Consider then a A^FG-dimensional foreground model, that is, the foreground vector is modeled as 
/ = Fg where P is a fixed A^obs x A^fg full column rank matrix and g is a. Nfg x 1 random vector. Then, just as in standard 
(one-dimensional) ILC, one may estimate g as g = W^^^x where W^^^ is an A^fg x A^obs matrix designed such that g 
has minimum power under the constraint of offering unit gain to the foreground. In other words, W^^^ is the minimizer of 
^|■^4^ILC^|2 ^j^^jgj. constraint W^^^F = In^g- The foreground emission is then reconstructed as f — Fg = FW^^^x. 
Using the Lagrange multiplier technique, it is readily found that 

/ = PfX where P/ = F{F' R'^ F)'^ F' R-\ 

Matrix P/ is recognized as an oblique projection onto Span(P) along the A^obs — A^fg dimensional subspace A//: 

J^f^{z\ F'R-^z = 0} . 

Hence, multi-dimensional ILC appears as a direct generalization of the one-dimensional case. 

It is interesting to determine the conditions such that P f — P2 because then the foreground estimation procedure by the 
CMB-ILC subtraction could be understood as a minimum variance procedure (the A^fg -dimensional ILC) and the minimum 
variance result would be obtained without even knowing matrix F. This is examined in the next subsection. 

A3 Subspace estimation 

Under which conditions do we have P2 = P/? First, we should have identical dimensions, i.e. Nfg = A^obs — 1- Let us then 
consider the noise-free model, cc = as + Fg where F is A^obs x (A^obs — 1) and Cov(gf) = G. The covariance matrix of the data 
is 

R = Cov(x) = aa' + FGF\ (Al) 

Given R and a, Eq. (Al) uniquely determines the range of F (see proof at the end of this section). Further, the corresponding 
value of P/ (which depends on F only through its range) is precisely equal to projection matrix P2. Therefore, we can conclude 
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that the CMB-ILC subtraction method is identical to the multi-dimensional foreground-ILC when the galactic subspace (the 
range of F) is determined from (Al). Note that these conclusions still hold when empirical values are used instead of ensemble 
averages. 



A3.1 Proof 

We show how the range space of F is determined from a and R — Cf^^aa^ + FGF^ when F has A^obs — 1 columns. 
Let W be an A^obs x A^obs whitening matrix i.e. satisfying WRW^ — I and denote whitened quantities with a tilde: 

X = Wx, a = Wa, F = WF. 

Consider the AT^bs x AT^bs matrix U = [a{Cf^^Y^^\ FG^^% Then UU' = Cf^^ ad' + FGF' = W{Cf^^ aa' + 
FGF')W' — WFiW' — I. Hence matrix U has ortho normal columns so that the range space of F is orthogonal to 
vector a. Therefore, the range space of F can be estimated as follows: take W any square root of ^ compute d — Wa, 
compute an A^obs x (A^obs — 1) matrix B whose columns are orthogonal to a so that range(S) is an estimate for range(l^), 
finally compute an estimate of range (i^) as range 



APPENDIX B: POWER SPECTRUM OF THE RECONSTRUCTION ERROR 

The temperature map at a given frequency channel of WMAP is given in harmonic space by 

Vim — OiSirn + + 5 

where i indexes the WMAP frequency bands, b} is the beam of the considered WMAP channel, and sim, flm ctnd are 
respectively the CMB, the foregrounds and the instrumental noise in the given frequency band. The CMB-ILC estimate at 
the W band resolution expands as 

S(.m — Oi S£rn + 0£m, 

where the residual noise of ILC (difference input-output) is 



Ei ( iW ni , i 



Here w^^c the needlet-ILC weights. Our foreground estimate is computed by subtracting the Wiener filtered CMB-ILC 
estimate at the beam of the considered frequency channel: 



The Wiener filter ae has been applied to reduce the residual noise of the CMB-ILC. It is thus given by 
where A^]^*^ = E [(^^m)^] • Therefore, the foregrounds estimate at the frequency channel i is 

m 1 

where the reconstruction error s}^ = /*^^ — h\fl^ expands as follows: 



The noise power spectrum Ei — E \^(e\^Y'^ at the frequency channel i (diagonal term of the covariance matrix of the error) 
is the sum of several contributions. We have: 

[(l-«.)&^]'^[5L] + (^) E[{5^mf]+E (n^^)' -2{l- a^)l^,(^^E[s^m6^ 



We neglect the last term (correlation between the CMB-ILC reconstruction error d^rn and the instrumental noise n^^) 
and we express the correlations between the CMB-ILC reconstruction error 5irn and the CMB sirn (loss of CMB power due 
to empirical correlations) by the following formula computed in Delabrouille et al. (2009) 



Sin 



= 4 X (l-^°b^) gCMB 



-0.02C> 
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where the number of observation frequencies is A^obs = 6 and the number of pixels is Np = 1024. One finds that 

2 



Ei = Nl + 



(1 - ai)^ + 26^ (1 - ai) [aiU^^ x O.i 



02 



^CMB 



where we note Nl = E [(n^m)^] , ^F^^ = E [{Simf] , and = E [sj^] . Introducing the following notation for the effect 

of both Wiener filtering and rebeaming 



web} 



we thus find 



Ei = Ni + rlNf-^ + 



(6^ (1 - a,))' + 26^ (1 - ae) bY n x 0, 



02 



^CMB 



REFERENCES 

Bouchet, F. R., & Gispert, R. 1999, New Astronomy, 4, 443 

Cardoso, J.-F., Le Jeune, M., Delabrouille, J., Betoule, M., & Patanchon, G. 2008a, IEEE Journal of Selected Topics in 

Signal Processing, vol. 2, issue 5, pp. 735-746, 2, 735 
Cardoso, J.-F., Martin, M., Delabrouille, J., Betoule, M., & Patanchon, G. 2008b, ArXiv e-prints 
Delabrouille, J., & Cardoso, J.-F. 2007, ArXiv Astrophysics e-prints 

Delabrouille, J., Cardoso, J.-F., Le Jeune, M., Betoule, M., Fay, G., & Guilloux, F. 2009, A&A, 493, 835 
Delabrouille, J., Cardoso, J.-F., & Patanchon, G. 2003, MNRAS, 346, 1089 
Efstathiou, G. 2004, MNRAS, 349, 603 

Fay, G., Guilloux, F., Betoule, M., Cardoso, J.-F., Delabrouille, J., & Le Jeune, M. 2008, Phys. Rev. D, 78, 083013 
Ghosh, T., Saha, R., Jain, P., & Souradeep, T. 2009, Phys. Rev. D, 79, 123011 

Gorski, K. M., Hivon, E., Banday, A. J., Wandelt, B. D., Hansen, F. K., Reinecke, M., & Bartelmann, M. 2005, ApJ, 622, 
759 

Guilloux, F., Fay, G., & Cardoso, J.-F. 2009, Appl. Comput. Harmon. Anal., 26, 143 

Herranz, D., Sanz, J. L., Hobson, M. P., Barreiro, R. B., Diego, J. M., Martmez-Gonzalez, E., & Lasenby, A. N. 2002, 
MNRAS, 336, 1057 

Hivon, E., Gorski, K. M., Netterfield, C. B., CriU, B. P., Prunet, S., k Hansen, F. 2002, ApJ, 567, 2 

Hobson, M. P., Jones, A. W., Lasenby, A. N., & Bouchet, F. R. 1998, MNRAS, 300, 1 

Komatsu, E., et al. 2009, ApJS, 180, 330 

Leach, S. M., et al. 2008, A&A, 491, 597 

Maino, D., et al. 2002, MNRAS, 334, 53 

Marinucci, D., et al. 2008, MNRAS, 383, 539 

Melin, J., Bartlett, J. G., k Delabrouille, J. 2006, A&A, 459, 341 

Melin, J., Bartlett, J. G., Delabrouille, J., Arnaud, M., Piffaretti, R., & Pratt, G. W. 2010, ArXiv e-prints 
Narcowich, F., Petrushev, P., & Ward, J. 2006, SIAM J. Math. Anal., 38, 574 
Patanchon, G., Cardoso, J.-F., Delabrouille, J., & Vielva, P. 2005, MNRAS, 364, 1185 
Remazeilles, M., Delabrouille, J., & Cardoso, J.-F. 2010, (in preparation) 
Saha, R., Jain, P., & Souradeep, T. 2006, ApJ, 645, L89 

Saha, R., Prunet, S., Jain, P., & Souradeep, T. 2008, Phys. Rev. D, 78, 023003 
Tegmark, M., & Efstathiou, G. 1996, MNRAS, 281, 1297 



14 T. Ghosh et. al 



K— band 




K-band - TOILC 




-0.30 ^^^^^^^^^K ^^^^^H 0.30 



K-band Foregrounds 




-0.30 ^^^^^A. ^^^^^H 0.30 



Figure Bl. K-band maps at K-band resolution 
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Figure B2. Ka-band maps at Ka-band resolution 
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Figure B3. Q-band maps at Q-band resolution 
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Figure B4. V-band maps at V-band resolution 
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Figure B5. W-band maps at W-band resolution 



